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Abstract: When plane waves diffract through fractal-patterned apertures, 
the resulting far-held prohles or diffractals also exhibit iterated, self-similar 
features. Here we show that this specihc architecture enables robust signal 
processing and spatial multiplexing: arbitrary parts of a diffractal contain 
sufficient information to recreate the entire original sparse signal. 
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1. Introduction 


Many natural systems exhibit fractal properties IQ ; in fact, scale invariance underlies many 
self-similar phenomena, from frost crystallization to animal colouration to stock-market pric¬ 
ing. In the realm of optics, the fractal anatomy of systems is widely associated with aggregated 
dielectric and metal colloids 0, crystals Q, and tissues 0. Fractal systems are also observed 
in nonlinear optics EH) and fractalized optical properties can even efficiently characterize or 
enhance the response of materials 00. Less well utilized are the features of diffractals, or 
the diffraction patterns of fractal signals |[9}|rT|. Diffractals feature in methods of encrypting 
data (ID as a versatile approach to double random-phase encoding GD However, to differen¬ 
tiate from previous work we exploit fractal architecture to improve transmission robustness. 

Here, we explore diffractals for their application in signal processing (^[10. We show 
that the free-space propagation of diffractal-signal architectures provides algorithmic value and 
spatial multiplexing properties; any arbitrary subsection of a diffractal contains sufficient infor¬ 
mation to recreate the original sparse signal that is transmitted with a specific fractal architec¬ 
ture. In a manner similar to compressive imaging (^[10 —where sparse signals reveal greater 
information via the diffraction through structures—here the fractal structuring within the signal 
sparseness prevents the loss of information. 

Like other applications of fractals in communications applications, the diffractal architecture 
exhibits trade-offs. Fractal antennas for the radio frequency and microwave regimes are known 
for being compact and versatile over wide spectral bands but are power intensive ; fractal 

encoding algorithms enable image compression with higher-resolution at the expense of greater 
algorithmic complexity p0| ; here, our research identifies that diffractal architectures prevent 
the loss of information but require greater signal preprocessing. This investigation extends our 
understanding of fractal structures in signal communications and may increase robustness and 
transmission rates of satellite, wireless, and interplanetary communication systems, i.e., net¬ 
works that support a large number of roaming receivers. 

The remainder of this report is organized as follows. First, we formalize the form of a trans¬ 
mitted fractal signal, show that its far-field diffraction pattern or diffractal is also a fractal, and 
illustrate the reciprocal nature of fractals with the Sierpinski carpet. Secondly, we demonstrate 
the robust retrieval of a signal from a diffractal; the original signal is reconstructed even when 
the majority of the diffractal signal is blocked. Finally, we discuss the future applications for 
diffractal spatial multiplexing in free-space communication systems and conclude. 


2. Theoretical Description 

2. 1 . Spatial Multiplexing of the Diffractal 

We generate the fractal transmittance pattern from any base matrix via recursive iterations 
where the matrix is resized and convolved with itself repeatedly pT| . The base matrix B{x^y) 







adopts a general form, 


( 1 ) 


Bi{.x,y) =Y^5{xr' ^-Xi,yr‘ ^-yj), 
j 

where the subscript denotes the iteration /, r > 1 is the relative scaling between iterations, and 
5{x — Xj,y — yj) is the Dirac delta function at x = xj and = yj- The fractal transmittance 
function T (x, y) is calculated recursively, 

T„{x,y) = Tn-i{x,y)*B„-i{x,y), (2) 

where the subscript denotes the order of the fractal or its expression at the n** iteration, To is 
the initial profile of I’s, and * denotes the convolution operator. 

Subsequently, the diffractal is the Fourier transform or the far-field of the transmittance func¬ 
tion f | [io| , 

n 

tn(kx-)ky) = tQ{kx-)ky)^^^Bi(kxtky) . (3) 

i=\ 

The Sierpinski carpet is one example of a fractal that is generated by this process and via the 
iterated substitution of a 3 x 3 base matrix of ones with removal of the center element: 
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( 4 ) 


The second substitution of Eq. represents the mathematical expression for the base ma¬ 
trix BQ{x^y). In the case of the Sierpinski-carpet base-matrix elements, r = 3, and xj and yj 
are the perimeter coordinates of a 3 x 3 9-unit block centered at the origin, and {xj^yj) G 
[(1,1), (1,0), (1, —1), (0,-1), (-1,-1), (—1,0), (—1,1), (0,1)]. Since each of the Dirac delta 
functions in 5o(-^5}^) yields a phase shift in the Fourier domain, i.e., ^{5{ax — Xi)} = 
^iKickxXi ^ 1 ^ 1 ^ scaled Fourier Transform at the iteration of the Sierpinski carpet [Eq. 
becomes: 


Bi{kx,ky) = {2/ry ^[cos(2;r3^ ^^;c)cos(2;r3^ ^^_y)+cos(2;r3 
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%)+cos(2;r3^-%)]. 


( 5 ) 


With each iteration, the spatial frequency components ky of the diffractal increase by a factor 
of 3 and spread the diffractal across a 3-times wider -range, which is evident in Eq.[^ the 
cutoff of fn scales in proportion with kx^kyoc . 

We calculate the Sierpinski carpet T recursively [Eq. : a second-order fractal is generated 


from the Kronecker tensor product of the base matrix 


fi 1 1 


with itself; a third-order fractal 
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is generated from the Kronecker tensor product of a second-order fractal and the same base 
matrix [see Code File 1] j^. Sierpinski carpets of ^ = 1,3, and 5 are shown in Fig.[T^) with 
the corresponding diffractals f [Fig.[^)]. 

As the fractal order increases, T exhibits smaller self-similar features at higher kx^ky\ the 
diffractal also exhibits a fractal architecture, as observed in other literature ©[I). Moreover, 
when n is large, an arbitrary subsection of t closely resembles the whole, and the self-similarity 
is already apparent with n = 5 [Fig.JT]:-d)]. The iterated, self-similar, wide-spatial-frequency 
features contained in the diffractal Tn [Eq. enable robust reconstruction of itself, which is 
described in the next subsection. 


2.2. Robust Reconstruction of a Blocked Diffractal Signal 

We now refer to the sparse matrix and transmitted signal T as the original and fractalized 
signal, OS and ES, respectively; the diffractal signal DS is T or the Eourier transform of ES; 









Fig. 1. a) Signal patterns, or fractalized signals (FS) of orders w = 1, 3, and 5. b) Corre¬ 
sponding Fourier transforms of signals on logarithm scale, or diffractal signals (DS). c) 
Reconstructed Fourier-transforms (RFS) from the 1% black-outlined subset, the blocked 
diffractal signal (BDS). d) Enlargement of a portion of the n = 5 diffractal, which illus¬ 
trates similar patterns at different length scales. 


BDS refers to an off-axis subsection of DS that is filtered; a reconstructed fractalized signal 
RFS refers to the inverse-Fourier transform of BDS; a regenerated version of the original signal 
ROS interpolates RFS in order to obtain OS. 

In Fig. [^), a subsection or BDS is outlined with a black square, which represents 1% of 
DS. The corresponding RFS from BDS are shown in Fig.[2c). For ^=1,3, RFS is ostensibly 
blank because BDS carries negligible power. In contrast, when n =5, RFS carries features that 
resemble FS. We refer to the capacity to reproduce FS with 1% of the off-axis and 0.1% of the 
transmitted power from the OS as robust reconstruction. 

Robust reconstruction—where RFS resembles FS—is possible even when the size of BDS 
is significantly reduced. A corresponding increase in the FS fractal order n will yield RFS that 
resembles FS when BDS is arbitrarily small. With the example of the Sierpinski carpet [Fig. [2, 
if the size of BDS is reduced to 0.1% of OS, then a comparable RFS is produced by increasing 
the fractal order of FS from ^ = 5 to ^ = 8. 

Robust reconstruction also refers to the capacity to regenerate OS from RFS from a simple 
threshold function [see Code File 2] | [2^ . The simple threshold function of the Sierpinski carpet 
divides RFS into a 3 x 3 array and measures the intensity in each of the 9 elements. Above a 
certain threshold value, the element is assigned a value of 1, and otherwise assigned a value of 
0. The 3x3 array that is processed with the threshold function, ROS, is theoretically identical 
to OS when OS is transmitted with the diffractal architecture. 

A partial explanation for the robust reconstruction is that increasing-order FS carry 
arbitrarily-high k^^ky and enable arbitrarily-^'ma// BDS to carry the information of FS or OS. If 
we strictly limit OS to binary or Dirac-delta functions, then FS has no kx^ky cut-off and as n ap¬ 
proaches infinity, fractalized features appear in FS without Sikx^ky cut-off. In fact, the amplitude 
of the additional kx^ky gained from each iteration scales inversely with which provides a 
detection limit only in practice; in theory, each iteration generates high-kx^ky copies of OS [Eq. 
1^ that are spatially distributed from the origin. 

Yet it is worth noting that the robust reconstruction is achieved because the diffractal archi- 
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Fig. 2. a) Light from a A = 632.8-™ wavelength laser is spatially filtered, expanded, and 
collimated and illuminates the full area of a 800x600 pixel spatial light modulator (SLM). 
The 4-f system is composed of 2 5-cm lenses placed after the SLM. The first lens Fourier 
transforms the fractal signal (FS) at the focal plane. An aperture is placed off-center at 
the focal plane and only transmits a portion of the diffractal (BDS). The second lens re¬ 
constructs the SLM image with the light that is transmitted through the aperture, b) The 
Sierpinski carpet hologram of order n= 5, CCD image in the focal plane of the hologram 
after the lens, and the reconstructed image. A circle denotes the area utilized to reconstruct 
the image. 


tecture also couples and ky in iterated products [Eq. [^ . As a result, RFS will resemble FS 
when either the BDS size or location changes. In a manner similar to spatial filtering, RFS will 
produce outlines of FS if n is not sufficiently large; however, unlike a high-pass spatial filter 
of a multi-scale random high-kx^ky pattern fTSl, a change in the location or the size of BDS 


will not distort the outline of RFS. Subsequently, the diffractal architecture provides superior 
performance over other algorithms that regenerate sparse data after filtering | T^p^ ; any of the 
reconstructed features in the ROS are strictly limited to the r x r-elements that are nonzero. 


3. Experimental Results 

We experimentally demonstrate the features of diffractals with a 4-/ optical arrangement 
where the 2-dimensional Fourier transform of a collimated fractal image, DS, lies in the focal 
plane of an imaging lens | [24| . The experimental setup that produces, filters and reconstructs 
FS is shown in Fig.|^) and experimentally reconstructed images, RFS, are shown in Fig.[^). 
An aperture of area approximately 0.8mm^ blocks the majority of DS. The placement of the 
aperture is shifted roughly 4mm horizontally and 5mm vertically from the central point, in the 
focal plane of the lens. The SLM image has a resolution of 800x600 pixels (16mmx 12mm). 

When we employ diffractal architectures of high orders {n >4), we observe the phenomenon 
that is shown numerically: the placement of the aperture in DS is irrelevant in order to recon¬ 
struct the original image. When the aperture is moved laterally in the focal plane, RFS maintains 
strong resemblance to FS. In fact, only the intensity of RFS diminishes as the aperture trans¬ 
lates further from the DS center; the outline remains fixed as the aperture moves. For smaller 
fractal orders, RFS resembles FS only when the aperture is placed within 0.8mm of the center, 
where the spatial frequency components are concentrated. The signal-to-noise of the experi- 





























































merit and CCD camera sensitivity limit effective reconstruction, while the highest order n is of 
FS is limited by the number of SLM pixels. 

4. Discussion of Applications 

Here we make the distinction that diffractals are specific fractal structures. Not all fractals 
enable robust signal communications and the recursive transformation employed to generate FS 
and DS in Eqs.[2 and [^differ fundamentally. For example, both of the Fourier Transform pairs 
FS and DS are fractals and carry iterated, self-similar features, but if we reverse their roles in our 
transmission system, the reconstruction will instead depend severely on the size and location 
of BDS. In fact, if in the example of the Sierpinski carpet, the center subsection becomes BDS, 
then the ROS will not resemble the OS, regardless of fractal order. The diffractal is unique 
from general fractals and self-similar scale-invariance alone is an insufficient precondition for 
our system of robust reconstruction and spatial multiplexing. 

It may seem contradictory that higher-order fractals lead to more robust signal transmis¬ 
sion since the finer structure of a higher-order fractal is itself harder to reconstruct. There are 
two perspectives of diffractals that explain the robust reconstruction. Firstly, the self-similar 
structures of higher-order fractals have a greater spatial frequency range and finer detail, and 
subsequently smaller BDS carry sufficient information to reconstruct OS. Secondly, the higher- 
order diffractals carry higher spatial-frequency components where the kx and ky components 
are coupled, and subsequently the location of the subsection in DS is unimportant. Here we 
have shown that BDS of arbitrary size and location carry sufficient information to reconstruct 
OS but in practice, there exists clear limitations for the robust reconstruction even in the limit 
of infinite-order FS. 

There exists a trade-off with the diffractal architecture between robust reconstruction and 
high bit rate; a greater bit-rate is achieved with a larger base matrix, which can limit the max¬ 
imum fractal order that is transmitted. For example, with the 3 x 3 or 9-element OS, there are 
512 possible spatial bits, three of which are illustrated in Fig. [^-c). A 4 x 4 base matrix re¬ 
quires more pixels than a 3 x 3 to achieve the same fractal order, n. With a limited SLM 
pixel resolution, there is a choice between the generation of higher-order fractals and the utility 
of higher spatial bits. 

If the trade-off between bit-rate and robust reconstruction are mitigated, then the diffractal 
architecture could support a large number of roaming receivers with only one transmitter, such 
as wireless or satellite networks shown in Fig.[^). The self-similar properties of the diffractal 
architecture and their corresponding far-field pattern provide a method to reach a large number 
of receivers, possibly moving, without signal degradation. The processing times required in the 
calculation of FS from OS are not trivial and scale with Secondly, the refresh rates of a 
spatial light modulator or similar adaptive-optics device present constraints on the maximum 
achieved bit rate, which requires further consideration. 

5. Conclusion 

We have shown that the diffractal architecture provides the beneficial features of spatial mul¬ 
tiplexing and robust reconstruction. We have centered our demonstrations with the Sierpinski 
carpet, a familiar fractal, although our results could have been demonstrated with 511 other 
patterns similarly imprinted with the diffractal architecture. Data that is transmitted with the 
diffractal architecture is highly robust to intermediate-obstacle signal blocks and diffractal sub¬ 
sections of arbitrary size and location carry sufficient information to regenerate the original 
signal without distorting outlines of its pattern. Our research illuminates potential applications 
in data transmission systems when one transmitter sends data to a large number of moving 
receivers or through noisy media. 



Fig. 3. Three examples of 512 9-bit spatial patterns, associated with base matrices a) 
[ 1 0 1 ] ’ b) [oil]’ [ 1 1 1 ] ’ fractal signals FS are shown with their correspond¬ 

ing experimentally-reconstructed fractal signals RFS from the experimental setup in Fig. 
[^). The lower-right inset shows the reconstructed original signal ROS. d) Example appli¬ 
cation: transmitted fractal signal FS is received at a far-held distance as a diffractal signal 
DS, where a roaming set of receivers, with only a diffractal subsection BDS, reconstructs 
the original signal OS. 
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